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Clarifying Forecasts of Dark Energy Constraints from Baryon Acoustic Oscillations 
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The measurement of baryon acoustic oscillations (BAO) from a galaxy redshift survey provides 
one of the most promising methods for probing dark energy. In this paper, we clarify the assumptions 
that go into the forecasts of dark energy constraints from BAO. We show that assuming a constant 
nPo.2 /G^ {z) (where P0.2 is the real space galaxy power spectrum at A: = 0.2ft/Mpc and redshift 
z) gives a good approximation of the observed galaxy number density expected from a realistic 
flux-limited galaxy redshift survey. We find that assuming nPo.2/G^(z) = 10 gives very similar dark 
energy constraints to assuming nPo.2 = 3, but the latter corresponds to a galaxy number density 
larger by ~ 70% at z = 2. We show how the Figure-of-Merit (FoM) for constraining dark energy 
depends on the assumed galaxy number density, redshift accuracy, redshift range, survey area, and 
the systematic errors due to calibration and uncertainties in the theory of nonlinear evolution and 
galaxy biasing. We find that an additive systematic noise of up to 0.4 — 0.5% per Az = 0.1 redshift 
slice does not lead to significant decrease in the BAO FoM. 
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I. INTRODUCTION 

It has been more than a decade since the discovery of 
the cosmic acceleration [l|, Illuminating the nature 
of dark energy, the unknown cause of the observed cos- 
mic acceleration, has become a prominent task for the 
astronomy and physics communities. 

Dark energy could be an unknown energy component 
3, or a modification of general relativity ^ . These two 
classes of models can be differentiated if both the cos- 
mic expansion history and growth history of cosmic large 
scale structure are accurately and precisely measured [5| . 
Current observational data are consistent with dark en- 
ergy being a cosmoiogicai constant, but other explana- 
tions are still allowed (see, for example, Q)- See [6| for 
recent reviews on dark energy research. 

One of the most promising methods for probing dark 
energy is to use the baryon acoustic oscillations (BAO) 
in the observed 3-D gala xy distribution as a cosmoiogicai 
standard ruler [sl. [ol. Il0| . At the last scattering of cosmic 
microwave background (CMB) photons, the acoustic os- 
cillations in the photon-baryon fluid became frozen, and 
imprinted their signatures on both the CMB (the acoustic 
peaks in the CMB angular power spectrum) and the mat- 
ter distribution (the BAO in the galaxy power spectrum). 
Because baryons comprise only a small fraction of matter, 
and the matter power spectrum has evolved significantly 
since last scattering of photons, BAO are much smaller in 
amplitude than the CMB acoustic peaks, and are washed 
out on small scales. BAO in the observed galaxy power 
spectrum have the characteristic scale determined by the 
comoving sound horizon at recombination, which is pre- 
cisely measured by the CMB anisotropy data [12]. In 
principle, the BAO scale can be extracted from data both 
in the transverse direction (s^), and along the line-of- 



sight (s||). Comparing the observed BAO scales with 
the expected values gives the angular diameter distance 
Da[z) = r{z)/(l + z) (where r{z) is the comoving dis- 
tance) in the transverse direction, and the Hubble pa- 
rameter H{z) in the radial direction: 



Da{z) 
S|| oc sH{z). 



(1) 



Seo & Eiscnstcin (2007) pJL] (henceforth SE07) pro- 
vided simple fitting formulae for forecasting the accura- 
cies of s/Da{z) and sH{z) from future galaxy redshift 
surveys. In this paper, we will clarify the assumptions 
made in BAO forecasts, and investigate their implica- 
tions. We will discuss the method in Sec. 2, present re- 
sults on BAO forecasts in Sec. 3, and summarize in Sec. 4. 



II. THE METHOD 

The simplest and most widely used method to forecast 
constraints from future observations is to use the Fisher 
matrix formalism. The Fisher information matrix of a 
given set of parameters, s, approximately quantifies the 
amount of information on s that we expect to get from 
our future data. The Fisher matrix can be written as 
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(2) 



where L is the likelihood function, the expected prob- 
ability distribution of the observables given parameters 
s. The Cramer- Rao inequality states that no unbiased 
method can measure the i-th parameter with standard 
deviation less than Xj yJFa if other parameters are known. 
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and less than \/ (F~^),;» if other parameters are estimated 
from the data as well |13l |. 

In the limit where the length scale corresponding to 
the survey volume is much larger than the scale of any 
features in the galaxy power spectrum Pg(k), we can as- 
sume that the likelihood function for the band powers 
of a galaxy redshift survey is Gaussian fl^. Then the 
Fisher matrix for estimating parameters from a galaxy 
redshift survey can be approximated as 



dpi 



^e//(k) 



2(27r)3 



(3) 



where pi are the parameters to be estimated from data, 
and the derivatives are evaluated at parameter values of 
the fiducial model. The effective volume of the survey 



n(r)Pg(fc,^) 



n(r)Pg(fc,/i) + l 



nPg{k,yi) 



nPg{k,fj.) + 1 



v.. 



survey ; 



(4) 



where the comoving number density n is assumed to only 
depend on the redshift for simplicity; /i = k • r/k, with f 
denoting the unit vector along the line of sight; k is the 
wavevector with |k| = k. Note that the Fisher matrix 
is the inverse of the covariance matrix of the parameters 
Pi if the Pi are Gaussian distributed. Eq.Q propagates 
the measurement errors in In Pg (k) (which are propor- 
tional to [Ve//(k)]~^/^) into measurement errors for the 
parameters pi. Note that Eq.([3]) can be rewritten as 



dPg{k,fl) dPg(k,fi) 



Pg{k,f,) 



dpt 

2TrPdk 
2 (27r)3 ' 



dpj 



where /x = k • r. 



A. The "wiggles only" fitting formulae 

SE07 [ll| provided simple fitting formulae for forecast- 
ing the accuracies of s/Da{z) and sH{z) from future 
galaxy redshift surveys. Essentially, they approximated 
Eq.® with 



Fij — ^survey 



' dfi dPi,ik,p\z)dPh{k,p\z) 

-1 Jk„i„ 



1 



P'™(fc,^|z) +n- 



dpi 
2TTpdk 
2(27r)3 ' 



dpj 



(6) 



where Pi,{k,fi\z) is the power spectrum that contains 
baryonic features. The linear galaxy power spectrum 



P^»(fc,Mk) 



p^:«(fciz)p(M) 



p^;;ik\z) = bizf 



G{z) 
G(0) 



(7) 

Plir{k\z = 0) (8) 



where P'™(fc|z) is the linear galaxy power spectrum in 



real space, h is the bias factor, G{z) is the growth factor, 
and P^^{k\z = 0) is the present day linear matter power 
spectrum. R{p-) is the linear redshift distortion factor 
given by 



PV) = (l + /3/^')' 



(9) 



The power spectrum that contains baryonic features, 
Pb{k,p,), is given by 



Pb{k,p\z) = V8^AoP';^{ko.2,t^\z) 



sin(a:) 



• exp 



\1.4 ^nl 



where we have define 



-1 



fco.2 = 0.2/iMpc 
= k • f = fc/i 



1/2 



ki_ = yjk^ = ky/l-p'^ 
The nonlinear damping scale 



(10) 

(11) 

(12) 
(13) 
(14) 



Si 



S±(l + /g 
12.4/1 



= 8.355/1 



•0.758 

G{z 



G{z) 



^--(S) "-c;(o) 



Pnl 



G(0) 



Pnl, 



(15) 



(5) where /, 



d\nG{z) / d\na denotes the growth rate of 
matter density fluctuations. The parameter pnl indi- 
cates the remaining level of nonlinearity in the data; 
with PNL — 0.5 (50% nonhnearity) as the best case, and 
Pnl = 1 (100% nonhnearity) as the worst case [ll|- For a 
fiducial model based on WMAP3 results [ll| (f^^ = 0.24, 
h = 0.73, 17a = 0.76, Vtk = 0, Vlbh^ = 0.0223, r = 0.09, 
ris , \ 

Silk damping scale Es — 8.38 h~^Mpc. 
Defining 



0.95, T/S = 0), Ao = 0.5817, P0.2 = 2710 cr|g, the 



Pi 

P2 



HDa/s), 



In s 



In sii = ln{sH) 



(16) 
(17) 



substituting Eq. (jTU]) into Eq. ^ , and making the approx- 
imation of cos^ X ^ 1/2, we find 



Fij — ^survey Aq 



d/^/i(M)/j(/^) 



dfc e 



^^"(fck = o) 1 

P}r{ko.2\z = 0) nP^™(fco.2,Mk)e~^-''^' 
exp [-2(fcS,)i-4-fc2sy ^ 



H -2 



(18) 
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where Pg™(fco.2, mI^) is given by Eq.d?]) with k = ko.2- 
Note that wc have added the damping factor, e^*^^'*^"'^, 
due to redshift uncertainties, with 



Or 

dz 



(19) 



where r is the comoving distance. For a flat universe, 
r — dz' / H{z'). The functions are given by 



/i(/i) = (91na;/91npi = /i^ — 1 
/2(m) = d\nx/d\TLp2= 

B. Full power spectrum calculation 



(20) 
(21) 



In order to compare the "wiggles only" method and 
the full power spectrum method for BAO forecast, we 
must include the nonlinear effects in the same way in 
both methods. In the full power spectrum method, the 
observed power spectrum is reconstructed using a par- 
ticular reference cosmology, including the effects of bias 
and redshift-space distortions 0: 



Pobs{kj_'^ , ku — 



[Da{z) 



[DA{z)rH{zY-f 



G{z) 



where 



and — 



G(0) 



Da{z) 
Da{zY^^'' 



PmatteriMz ^ 0) + Pshot, (22) 



fcl 



H{z) 



ref 



H{z) 



(23) 



kyk^ = fcjj/(fcj^ + A:j|). The values in the 
reference cosmology are denoted by the subscript "ref" , 
while those in the true cosmology have no subscript. 

Following Seo & Eisenstein (2007) we include non- 
linear effects in the full power spectrum calculation by 
modifying the derivatives of P(k) with respect to the pa- 
rameters Pi, i.e.. 



dPg{k,tJi\z) _ dP^"'ik,fi\z) 



dpi 

Eq.® becomes 



dpi 



exp {~k^J:lJ2) . (24) 



nP'"'{k,p) 



ik,^i) + l 



91nPj™(A;,At) ainPj*"(fc,/i) 
dpi dp J 

_j,2j^2 2'!rk^dk 
2(27r)3 ■ 



confidence contour for the parameters {wq, Wa), assuming 
a dark energy equation of state wxiz) — wq + {I — a)wa 
(TtI . It is most convenient to define a relative generalized 
FoM M 



FoM.. 



1 



VdetCov(/i,/2,/3,...)' 



(26) 



where {fi} are the chosen set of dark energy parameters. 
This definition has the advantage of being easy to calcu- 
late for either real or simulated data. When applied to 
(wo, Wa), we find 



FoM^ = 



1 



1 



^det Cov(wo, Wa) 

(27) 

which differs by a factor of 6.177r from the DETF defi- 
nition. This is the FoM that has been widely used, and 
tabulated in the DETF report. 

The most sensible FoM requires choice of dark energy 
parameters that are least correlated [l^. A good choice 
is to use (wo,Wo.5) from 



wx{a) = 3wq,5 - 2wq + 3{wo- W0.5) a 



(28) 



X{z) = (1 + zf{l-2wo+3wo.,) 



9{wo - W0.5) 



1 + z 



where wq.s is the value of wx at z = 0.5. The correla- 
tion of (wojWq.s) is much smaller than that of (wo^Wa)- 
For real data analyzed using Marcov Chain Monte 
Carlo (MCMC), the covariance matrix of (wojWo.s) and 
(wo, Wa) cannot be transformed, since choosing (wo, W0.5) 
and choosing (wq, Wa) as the base parameters correspond 
to different priors, if uniform priors are assumed for the 
base parameters. For Fisher matrix forecast, one can 
simply transform the covariance matrix of (wQjWa) into 
that of (wojWo.s) by using 



Wo. 5 = Wo+ Wa/i- 



D. Choice of redshift bins 



(29) 



There has been some confusion in the literature about 
the optimal choice of redshift bins in constraining dark 
energy. This choice should depend on the observational 
method considered. For a galaxy redshift survey, the 
observables are s/Da and sH (length scales extracted 
from data analysis). Since these scales are assumed to be 
constant in each redshift slice, the redshift slices should 
be chosen such that the variation of 1/Da{z) and H{z) 
(25) in each redshift slice remain roughly constant with z. 

For a flat universe, the variations of H{z) and 1/ Da{z) 
with z are 



C. Figure of merit and dark energy parameters 



diH/Ho) _ 3nr„{l + zY + nxX'{z) 



The DETF defined the dark energy figure of merit 
(FoM) to be the inverse of the area enclosed by the 95% 



dz 
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dz 
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r{z) 
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2E{z) 

l + z 
H{z)r{z)_ 



(30) 
(31) 
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FIG. 1: The variations of H(z) and 1/Da{z) with z. 

where X{z) = px{z) / px{0) is the dark energy density 
function. X{z) = 1 for a cosmological constant. Fig[T] 
shows the variations of H{z) and 1/Da{z) with z for a 
fiducial flat ACDM model. The amplitude of variation 
in H{z) increases slowly with z. The amplitude of varia- 
tion in 1/ Da{z) decreases with z, and stablize for z > 1. 
Thus choosing a constant Az Jar redshift slices is the op- 
timal choice. Choosing constant Az avoids significantly 
degrading the approximation of H{z) being constant in 
each redshift slice as z increases, while the approximation 
of Da{z) being constant becomes increasingly better as 
z increases. This makes sense since H{z) measurements 
carry more weight in constraining dark energy than the 
Da(,z) measurements. 

III. RESULTS 

We find that the fitting formulae, Eg. pS)) . give errors 
of hi{s/DA) and In(siJ) that match those from the full 
Fisher matrix, Eq.(Iin]), to better than - 10-20%. We 
have applied the fitting formulae, Eq. fTS)) . to investigate 
the dependence of the forecast of dark energy constraints 
from BAO on various survey parameters. We adopt the 
Dark Energy Task Force (DETF) fiducial model, with 
= 0.146, VLbW^ = 0.024, h = 0.725, Ofc = w; = -1, 
and ris = 1.0 [Hi. We assume cr8,m(^ = 0) = 0.8 [13], 
and (T8,.g{z = 0) = 1. We compute the matter transfer 
function using CMBFAST 4.5.1 [H. Table 1 lists the 
parameter values in the DETF model and the "other" 



model for comparison (all are present day values). We 
take krnin = 10^'*/iMpc~^, and k^ax — 0.5/iMpc~^ [ll|- 



A. Dependence on galaxy number density 

For a given BAO survey, the galaxy number density 
n{z) and bias function b{z) should be modeled using 
available data and supplemented by cosmological N-body 
simulations that include galaxies [27| . Since n{z) and 
h{z) depend on survey parameters such as the flux limit 
and the target selection method, a more generic galaxy 
number density given by assuming nPo. 2 = nPg{kQ,2\z) — 
constant is often used in BAO forecasts, where Pg{kQ,2\z) 
is the real space power spectrum of galaxies at fc = 
0.2ft,Mpc~^ and redshift z. Fig [2] shows the number 
of galaxies per Az = 0.1 slice per unit volume and 
per (deg)^ respectively, multiphed by 6^(z). Clearly, as- 
suming a constant nPo.2 corresponds to assuming that 
n(z)b'^{z) increases faster than a linear function in z, 
most likely faster than the increase in the galaxy bias 
function h{z) [13, H^, thus implying an assumed galaxy 
number density n{z) that increases with redshift. 

A typical flux-limited galaxy redshift survey has a 
galaxy number density distribution that peaks at some 
intermediate redshift (as the volume per redshift slice 
increases with z), and decreases toward the high z end 
of the survey (as the increase in the number of galaxies 
fainter than the flux limit dominates over the increase in 
volume per redshift slice). For a galaxy redshift survey 
with limited resources, photometric pre-selection could 
be used to alter n{z) to cut out low redshift galaxies, 
and validate the assumption of constant riPo.2- However, 
this would correspond to nPo.2 < 1- It is challenging to 
achieve nPo.2 > 1 at z = 2 even for the most optimistic 
space-based galaxy redshift survey [HI, HI]. Thus it is 
misleading to use nPo.2 = 3 to forecast the dark energy 
constraints from a space-based survey. 

We propose the following alternative assumption about 
the observed galaxy number density: 

nPQ,2/G'^{z) ^nb^{z)P„^{kn,2\z = 0) = constant. (33) 

FigO shows the expected rtPo.2 assuming constant 
"■^b.2/G^(z) in the DETF model. Assuming a constant 
nPo.2/G'^(z) allows us to avoid assuming an unrealisti- 
cally high riPo.2 at the high redshift end, while at the 
same time retain the increasingly higher nPo.2 easily 
achievable at intermediate and lower redshifts that play 
a key role in tightening the dark energy constraints. 

In forecasting dark energy FoM from BAO, it is im- 
portant to specify whether a constant nPo.2 or con- 
stant 7iPo.2/G^(z) is assumed. FigH] shows the frac- 
tional errors of {sH) and {s/Da) per Az = 0.1 slice ex- 
pected from a 28,000 (deg)^ galaxy redshift survey with 
f^z/i^ + z) ~ 0.001, using the "wiggles only" fitting for- 
mulae of Eq. (fT5)l , for different assumptions about the ob- 
served galaxy number density. Note that for the "other" 
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TABLE I: Two fiducial cosmological models 



name 
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Us (78, m 




DETF 
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0-006 
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FIG. 2: The comoving number density of galaxies per unit co- 
moving volume (upper panel) and per (deg)^ (lower panel). The 
results for two different fiducial cosmological models are shown, the 
"other" model (solid), and the DETF model (dotted). 



model, our results are in excellent agreement with the re- 
sults of SE07 [ll| for nPo.2 = 3 (thin lines with squares, 
compare with the dashed lines in Fig. 3 of SE07), and for 
the cosmic variance case, i.e., zero shot noise and zero 
nonlinearity (thin lines with triangles, compare with the 
sohd lines in Fig.3 of SE07). 

Figini shows the 68.3% joint confidence contours for 
{wo,Wa) and (^0,^0.5) (see Sec III CI for the definition 




0,5 1 1.6 2 

z 



FIG. 3: The expected nPo.2 assuming constant nPo,2/G'^{z) in 
the DETF model. 



of the parameters) for the three different galaxy num- 
ber densities shown in Figl2]and FigHl The DETF fidu- 
cial model is assumed. It is interesting to note that as- 
suming TiPo. 2/^^(2:) = 10 gives very similar dark energy 
constraints to assuming rtPo.2 — 3, but the latter corre- 
sponds to a galaxy number density larger by ~ 70% at 
z = 2 (see Fig^. 

Figini shows the relative dark energy FoM, FoM^ (de- 
fined in Sec lII Cp . for a galaxy redshift survey covering 
28,000 (deg)^ and 0.3 < z < 2.1, as a function of galaxy 
number density. The dashed lines in FiglS] show the re- 
suhs of using the full P{k) method (see Eq.[l5]). The full 
P{k) method boosts the FoM by a factor of ^3-4, with 
more gain for higher galaxy number densities. 



B. Dependence on other survey parameters 

We now study the dependence of the dark energy 
FoM on the other survey parameters, redshift accuracy, 
redshift range, and survey area. We show all our re- 
sults for two representative galaxy number densities, 
nPo.2/G'^iz) = 3 and 10. 
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A[sH] 
sH „ 





• true model 




- BAO [wiggles only, 28.000 (deg)^ cr,/(l+z) = 
Nonlinearity=50%] 
BAO + Planck 



FIG. 4: The fractional errors of (sH) and (s/Da) per Az = 0.1 
slice expected from a 28,000 (deg)^ galaxy redshift survey with 
0-2/(1 + 2) = 0.001, using our "wiggles only" modified formula 
of Eq.l ll8l l. for different assumptions about the observed galaxy 
number density The results for two different fiducial cosmological 
models are shown, the "other" model (solid), and the DETF model 
(dotted). 
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0.2 


nP 

nP 
nP 


,/G'(z).3 
,/G>(z).l0 

2=3 









-0.2 


IDJao 






BAO [wiggles only, 28.000 (deg)=. 

Nonlinearity=50%] 

BAO + Planck 

!...!. 


»y(i+z)..ooi, 

1 . . 



Fig[7| shows the FoM^ for a galaxy redshift survey cov- 
ering 28,000 (deg)^ and 0.3 < z < 2.1, as a function 
of the redshift accuracy. The FoM decreases quickly as 
az/{l + z) increases beyond 0.001. Increasing the redshift 
accuracy to CTz/ + < 0.001 does not have a significant 
impact. 



FIG. 5: The 68.3% joint confidence contours for (u)o,tUa) and 
{wo, wo. 5) for the three different galaxy number densities shown in 
Fig[2]and Fig|4] The DETF fiducial model is assumed. 



Fig [8] shows FoMr for a galaxy redshift survey covering 
28,000 (deg)2 and (tJ(1 + z) = .001, as a function of 
Zmax for Zmin = 0.5 (upper panel), and as a function of 
Zrmn for z^ax = 2 (lowcr panel). Decreasing Zmax and 
increasing Zmin both significantly decrease the FoM for 
BAO only. 

Fig[5] shows the FoM for a galaxy redshift survey cov- 
ering 0.3 < z < 2.1 with az/ {1 + z) — .001, as a function 
of survey area. The BAO only constraints scale linearly 
with survey area. This is as expected, since no priors 
are added, and the constraints on dark energy parame- 
ters scale with l/-y/area. The more external priors are 
added, the more slowly the FoM increases with survey 
area. Thus the FoM of (BAO-|-Planck) increases more 
slowly with area than that of BAO only, and the FoM of 
(BAO-|-Planck+StageII) increases with area even more 
slowly. 



C. Dependence on systematic uncertainties 

There are two types of systematic uncertainties for 
BAO. One type of uncertainty is due to the calibration 
of the sound horizon by CMB data. The other type of 
uncertainty is due to the intrinsic limitations of the BAO 
method, for example, the uncertainties in the theory of 
nonlinear evolution and galaxy biasing [2^ . 

WMAP five year observations gives As/s ^ 1.3% (2^ . 
We expect that Planck wiU give As/s ^ 0.2% Note 
that s appears as an overall scale parameter in the observ- 
ables [s/Da) and [sH), and is assumed to be statistically 
independent from {s/Da) and {sH). Thus the uncer- 
tainty in s has no effect on the dark energy FoM for the 
BAO wiggles only method. However, when other data are 
combined with the BAO (wiggles only) data, the overall 
dark energy FoM is decreased due to the uncertainty in 
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BAO: aa.OOO (deg)«, 0.3<z<2.1, +z) = .001 

FoMj, for wiggles only 

[FoM, tor full P(k)]/3 




nP„VG=(z) 

FIG. 6: The relative dark energy FoM, FoMr, for a galaxy redshift 
survey covering 28,000 (deg)2 and 0.3 < z < 2.1, as a function 
of galaxy number density. Both BAO wiggles only method (solid 
lines) and the full P{k) method (dashed lines) are shown. The 
DETF fiducial model is assumed. 



s. If the parameter set used is {wo,Wa,^x,^ki^mh^) 
for the BAO dark energy covariance matrix, then As 7^ 
only modifies the diagonal matrix for flmh^ by adding a 
term 



^DE.55 



(34) 



The Planck prior of As/s ~ 0.2% has a very small effect 
on the FoM. 

Our current modeling of intrinsic BAO systematic ef- 
fects are at the ~l-2% level for realistic N-body simula- 
tions that include galaxies (not just dark matter or dark 
matter haloes) [131 • Correcting the bias in the estimated 
BAO scale due to systematic effects will likely lead to 
an increase in the uncertainty of the derived BAO scale 
[Is*! . In order to make realistic assessment of dark energy 
constraints from BAO, one should allow for some level of 
remaining system uncertainty in each redshift bin. 

We show the effect of systematic uncertainty in two dif- 
ferent ways. First, we show the effect of nonlinear effects 
explicitly, using the parametrization shown in Eq. ()15|) . 
with pnl indicating the level of remaining nonlinear- 
ity. Secondly, we will show the effect of additive sys- 
tematic noise in each redshift slice due to the incom- 
plete removal or imperfect modeling of nonlinear effects 
or scale-dependent bias. 

Fig [6] shows the relative dark energy FoM for three lev- 
els of nonlinearity: 50%, 75%, and 100%. Clearly, dark 
energy constraints from BAO are extremely sensitive to 
the level of nonlinearity assumed. FigfTUlshows the FoM^ 
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FIG. 7: The relative dark energy FoM, FoMr, for BAO wiggles 
only from a galaxy redshift survey covering 28,000 (deg)^ and 0.3 < 
2 < 2.1, as a function of the redshift accuracy. The DETF fiducial 
model is assumed. 



as a function of the level of nonlinearity for two repre- 
sentative galaxy number distributions. 

FigdU shows the FoM^, for a galaxy redshift survey 
covering 28,000 (deg)^ and 0.3 < z < 2.1, as a func- 
tion of the level of additive systematic errors in each red- 
shift bin. The pessimistic case considered by the DETF, 
aly, = 0.01 X gives 0.0224 for Az = 0.1, be- 

yond the range of Fig llll An additive systematic error 
of 0.5% per Az = 0.1 redshift slice would decrease the 
BAO FoM by 23% for nPo.2/G^{z) = 3, and 34.5% for 
nPQ.2l G'^{z) = 10. Note however, the larger galaxy num- 
ber density should allow better modeling of the system- 
atic effects, thus the systematic errors should be lower for 
riPo.2/G^(z) = 10. An additive systematic error of 0.35% 
per Az = 0.1 redshift slice would decrease the BAO FoM 
by 21% for nPQ,2/G'^{z) = 10. 
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FIG. 8: The relative dark energy FoM, FoMr, for BAO wiggles 
only from a galaxy redshift survey covering 28,000 (deg)^ and 
(Tz/(1 + z) = .001, as a function of Zmax for Zmin = 0.5 (upper 
panel), and as a function of z^nin for Zmax = 2 (lower panel). The 
DETF fiducial model is assumed. 



IV. SUMMARY 

We have studied the various assumptions that go into 
the BAO forecasts, using the fitting formulae for BAO 
("wiggles only") We have shown that assuming 

nPo.2/G'^{z) = constant gives a more realistic approxima- 
tion of the observed galaxy number density than nPo.2 = 
constant. We find that assuming 7iPo.2/G^(z) ~ 10 
gives very similar dark energy constraints to assuming 
nPo.2 = 3, but the latter corresponds to a galaxy num- 
ber density larger by ~ 70% at z = 2 (see FigIS]). 

Assuming that nPo.2/G^{z) = constant, we have shown 
how the FoM for constraining dark energy depends on 
the assumed galaxy number density, redshift accuracy, 
redshift range, survey area, and the systematic errors 
due to uncertainties in the theory of nonlinear evolution 
and galaxy biasing. 

We find that (Jz/(1 -I- = 0.001 is an optimal redshift 
accuracy for a galaxy redshift survey, with the BAO FoM 



.' ' 1 ' 
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BAO: wiggles only, 28,000 (deg)', ci,/(l+z)=.001 
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FIG. 10: The relative dark energy FoM, FoMr, for BAO wiggles 
only from a galaxy redshift survey covering 28,000 (deg)^ and 0.3 < 
2 < 2.1, as a function of the level of nonlinearity. The DETF 
fiducial model is assumed. 
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0-^, is added in quadrature to iln(3/D,) 
and Aln(sH) per Az=0.1 slice 



BAO: wiggles only, 28.000 (deg)=, Nonlinearity=5D%, is/s = 0.2% 




(7^, is added in quadrature to Aln(s/D,) 
and Aln(sH) per Az=0.1 slice 



FIG. 11: The relative dark energy FoM, FoM,-, for BAO wiggles 
only from a galaxy redshift survey covering 28,000 (deg)^ and 0.3 < 
z < 2.1, as a function of the level of systematic errors in each 
redshift bin. The DETF fiducial model is assumed. 

decreasing sharply with increasing Ozj (1 + z) for (T^/(1 + 
z) > 0.001. Further, the BAO FoM is very sensitive to 
the redshift range of the survey (see Fig[5]), at both the 
low redshift and high redshift ends. The optimal galaxy 



redshift survey should measure the BAO using the same 
tracer over the entire redshift range in which dark energy 
is important, i.e., < z < 2, this would enable robust 
modeling of systematic effects, as well as enabling the 
strongest dark energy constraints from BAO only. 

The FoM of BAO is very sensitive to the level of non- 
linearity assumed. If the level of remaining nonlinearity 
is 80%, instead of the best case of 50%, the FoM of BAO 
is decreased by almost a factor of two (see Fig|TU|. The 
assumed sound horizon calibration error of As/s = 0.2% 
(expected from Planck) has no effect on BAO only con- 
straints, and a negligible effect on the FoM of BAO com- 
bined with other data. An additive systematic noise of 
up to ^ 0.4 — 0.5% per Az — 0.1 redshift does not lead 
to significant decrease in the BAO FoM (see Fig fTTj) . 

Finally, we note that future dark energy surveys are 
usually assessed by their performance when combined 
with Planck and Stage II priors [To^. When this is done, 
it should be clearly stated. Adding Planck priors to BAO 
typically boosts the FoM by about a factor of ^ 5 com- 
pared to BAO only. Adding both Planck and DETF 
Stage II priors boosts the FoM by about a factor of ~ 10 
compared to BAO only (see FigsiTllTTI). FigH FigfTOl 
and FigHT] show that the FoM of combining BAO with 
Planck is less sensitive to redshift errors and systematic 
errors. This is more so when Stage II priors are added. 
The BAO FoM (BAO alone or combined with other data) 
also depend the choice of the fiducial cosmological model. 

A transparent comparison of the forecasts from differ- 
ent BAO projects can only be achieved if all are required 
to choose (in the order of impact) (1) the same BAO ap- 
proximation method (for example, the wiggles only fit- 
ting formulae of Eq.[TH]); (2) the same priors from Planck 
and current or ongoing projects (for example, the DETF 
Stage II priors); (3) the same fiducial model (for exam- 
ple, the DETF model, or the five year WMAP bestfit 
model) . 
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